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Abstract. The basic model of the non-equilibrium low dimensional physics the so-called 
totally asymmetric exclusion process is related to the 'crystalline limit' (q — ► oo) of the 
C*~) ' SU q (2) quantum algebra. Using the quantum inverse scattering method we obtain the exact 

expression for the time-dependent stationary correlation function of the totally asymmetric 
simple exclusion process on a one dimensional lattice with the periodic boundary conditions. 
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1 Introduction 

The study of quantum exactly solvable models within the framework of the quantum inverse 
OO | scattering method (QISM) [TJ El [3] has led to an algebraic structures which are g-deformations 

^sO ■ of the universal enveloping algebras [H [5j. These deformed structures are usually called now 

'quantum groups' or 'quantum algebras'. Integrable models associated with the 'crystalline 
limit' (q — > oo) [6] of these algebras play an important role in enumerative combinatorics [8] 
' and non-equilibrium low dimensional physics [U [10] . 

<^ • The algebraic approach to the calculation of the correlation functions for the integrable 

models within the frames of the (QISM) was developed in the late eighties of the last century 
advantage of this approach is that the form of the final answers depends on the i?-matrix of the 
model and not on the particular realization of the operators intertwained by this matrix. It is 
also important that the operator averages are calculated over the eigenstates of the generating 
function of the integrals of motion. 

The correlation functions for the boson model associated with the 'crystalline limit' of the 
XX Z i?-matrix were calculated in [12\ [T3] where it was shown that they have the determinantal 
form. 

The totally asymmetric simple exclusion process (TASEP) is an integrable model that is 
connected with the 'crystalline limit' of the XX Z i?-matrix. It is one of the most studied models 
of the low dimensional non-equilibrium physics. The Bethe ansatz approach to the solution of 
the model was developed in the papers [HI HSJ HH H3 HHJ H9]. The relation of the model to 
the QISM has been studied in [201 ED E2]j the model with the twisted boundary conditions 
within the frames of this method was considered in [23], and the asymmetric exclusion process 
in [21]. The dynamical properties of TASEP are intensively studied (see, e.g., [251 ESI |27] and 



*This paper is a contribution to the Proceedings of the XVIIth International Colloquium on Integrable Sys- 
tems and Quantum Symmetries (June 19-22, 2008, Prague, Czech Republic). The full collection is available at 
|http:// ww.emis.de/journals/SIGMA/ISQS2008.html | 



2 



N.M. Bogoliubov 



the references in these papers). The KPZ scaling limit of TASEP was considered, in particular, 
in HUES ESI- 

In the present paper we shall consider the TASEP on a periodic lattice with the finite number 
of sites and with the finite number of particles. We shall represent the form-factors of the local 
operators in the determinantal form using the algebraic approach of the QISM. This result may 
be considered as the main result of the present paper. 

The paper is structured as follows. In Section [21 the definition of the model is given. The 
algebraic approach to the solution of the model is presented in Section [3j in this section the 
stationary states are expressed through the state vectors of the model and the scalar product 
of the state vectors is expressed in the determinantal form. In Section HI the calculation of the 
averages of the local projection operators over the state vectors of the model defined on a lattice 
with M sites is reduced to the calculation of the scalar product of the state vectors of the model 
defined on the M — l sites. In Section [5j the well known results on the Bethe ansatz solution are 
presented. The answer, based on the results of the section [6l for the time-dependent stationary 
correlation function is obtained. Finally, conclusive remarks are given in Section [71 

2 Totally asymmetric simple exclusion model 

The totally asymmetric exclusion process describes a system of N particles on a periodic ring 
with M sites labelled i = M, M — 1, ... ,2,1. The particles move randomly only in one direction. 
We shall choose this direction from right to left. The exclusion rule forbids to have more than 
one particle per site. 

This process can be conveniently represented using the spin description - in each lattice site 
the spin- up state corresponds to the empty site |0), the spin-down state corresponds to the 
occupied one |1). The configuration of iV particles (2N < M) located on the sites M > m\ > 
m2 > • • • > m N > 1 is given by the basis vector 

\m 1 ,m 2 , . . .,m N ) = a mi a m2 ■ ■■a~ N \n), 

where the generating state is the state with all spins up 

|n> = o^ 1 |o>,. 

The described process is defined by a non-Hermitian Hamiltonian [14j : 

3=1 k J 

Here o" ±,z are the Pauli matrices, the matrix with subindex j acts nontrivially only in the j-th 
spin space of the total space of states of the chain (C 2 )® M : Sj = I®---®I®s®I®---®I, 
and the periodic boundary conditions are assumed: a n = a n +M- 

In the matrix representation —H is the Markov matrix that controls the time evolution of 
the probability measure. 

The first term of the Hamiltonian describes the hoppings of particles 

o7+i°i"l >i+il 1 )i = l 1 )i+il°)i- 

The second term counts the number of the allowed jumps of the particles with hard-core repul- 
sion. 

In our paper we shall consider the time-dependent stationary correlation function [T2] 

Z^iSN^e-^SmlSN), (2) 
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where the projection operator 

Sk = \{l + o z k ) (3) 

has value 1 if there is no particle at site k and otherwise. In the steady state |SV) every spin 
configuration with iV spins down has the equal weight: 

\Sn) = \mi,...,m N ). 

M>mi>m2>->mpf>l 

This state is a ground state of Hamiltonian with the eigenvalue zero: 
H\S N ) = 0. 

The left eigenvector of H with the eigenvalue zero is 
(S N \= ^ (m 1 , . . . ,m N \. 

M>mi>m2>--->mjv>l 

The total number of configurations in the ground state 



Zn = {Sn\Sn) 



N\(M -N)\ 

One can express the right steady state in the following form [30| : 

\s N ) = p%\n), (4) 

where is the iV-th power of operator 

M 



'Pr = '^2 Sm ' ' ' Sk + ia k 



k=l 

Sj are projection operators ([3]) and are Pauli matrices. Similarly for the left steady state one 
has 

(s N \ = (n\pg (5) 

with the operator Vl equal to 

M 



V L = ^o-^s fc _i ■■■s 1 . 



k=l 



3 Solution of the model 

To the solution of the model we shall apply the quantum inverse scattering method [U El [3]. 
The L-operator of the considered model \21\ [23] is a 2 x 2 matrix with the entries acting on the 
space of states of an M-site spin-^ chain: 

L(n\u) = ( u] . ° T ^-i Sn 1 =ss n + (I - s)(ul - n _1 s n ) + a~a+ + cr+a', (6) 

where the parameter u G C. 
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The L-operator ([6]) satisfies the intertwining relation 

R(u, v) (L(n\u) <S> L(n\v)) = (L(n\v) <g> L(n\u)) R(u, v), 
in which R(u, v) is the 4x4 matrix 



R(u, v) 



( f(v,u) 

g(v,u) 1 

g(v,u) 

V f(v,u) J 



(7) 



with 



f(v,u) 



uv 



9 9 ' 

vr — v z 



g{v,u) 



,2' 



In the quantum group theory this i?-matrix is known as the 'crystal base' i?-matrix [6]. 
The monodromy matrix is the product of L-operators 



T(tt) = L(M\u)L(M - l\u) ■ ■ ■ L(l\u) 



A(u) B(u) 
C{u) D(u) 



(8) 



The commutation relations of the matrix elements of the monodromy matrix are given by the 
same i?-matrix JTl) 



R{u, v) (T(n) T(v)) = (T(v) T(u)) R(u, v). 

In the explicit form the most important of these relations for us are 

C(u)B(v) = g(u, v) {A(u)D(v) - A{v)D{u)} , 
A(u)B(v) = f(u, v)B{v)A{u) + g(v, u)B(u)A(v), 
D(u)B{v) = f(v,u)B(v)D(u) + g(u,v)B(u)D(v), 
[B(u),B(v)] = [C(u),C(v)]=0. 

The transfer matrix t(u) is the matrix trace of the monodromy matrix 

T (u) = u- M HT{u) = u~ M (A(u) + D(u)) . 



(9) 



(10) 



(11) 



The relation ([9]) means that [t(u),t(v)] = for arbitrary values of parameters u, v. 

The cyclic shift operator r in the total spin space (C 2 ) ® M is expressed through the transfer 
matrix: 



t = t(1) 



n 12 n 23 ■ • • n 



M-1M- 



(12) 



Here 



II m n — S m S n + (/ S m )(J Sn) + (T m (7 n + (T m CT n , 

is the permutation operator: IL mn a m = a n Yl mn . The shift operator shifts the site indices 

T G\T = <J n 

and possesses the property t m = I. 

In terms of the transfer matrix (jlip the Hamiltonian (JTJ) is given by 



(13) 
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The right state vector of the model is the vector generated by the multiple action of operators 
B(u) = u~( M ~^B(u) on the generating state \Q) = ®-^ 1 |0)j 



N 

\y(u l ,u 2 ,...,u N ))=l[B(u i )\n). (14) 

i=l 

The generating state is annihilated by the operator C{u) 

c(u)\n) = o, 

and it is an eigenvector of operators A{u) and D(u), 

A( u )\n) = a (u)\n), D(u)\n) = s(u)\n) (15) 

with the eigenvalues 

a(u) = u M , 6(u) = (u - u-y 1 . 
The left state vector is equal to 

N 

(*(«!, «2,...,«JV)| = (n\Y[C( Ui ), (16) 

i=l 

where C( Ui ) = u~^ M ~^C{u), and (Sl\B{u) = 0. 

From definitions ([6]) and (jSJ) one finds by a direct calculation that u M ~ l B{u) and u M ~ 1 C(u) 
are polynomials in u 2 of degree M — 1: 

u M ~ l B(u) = u^-^Vr + ■■■ + (-l) M -V^ M -i ■ ■ ■ -i, 
n M " 1 C( U ) = u 2 ^- 1 )^ + • • • + (-l^-V/ • • • S2 a+ 

This decomposition allows us to express the steady states Q and © through the state vec- 
tors (]14p and (I16|) . namely 



TV TV 



lim Ts(«i)|n), (5jv|= lim (n|TTC(«i). (17) 

{«}— »oo "~ {u}— »oo 



The scalar product of the state vectors (|14p and (|16p may be evaluated directly by means 
of the commutation relations (|10|) with the help of Laplace formula for the determinant of the 
sum of two matrices. It also may be obtained as a special limit of the general formula for the 
scalar products (IX. 6. 26) of [3]. It should be mentioned that for the i?-matrix (|7|) the auxiliary 
quantum fields (IX. 6. 24) commute and the answer for the scalar product has the determinantal 
form. It means that the scalar products of the state vectors of the integrable models associated 
with the i?-matrix (J7|) are expressed as the determinants of the matrices with the entries which 
depend on / and g elements of i?-matrix and on the eigenvalues a and 5 (|15p that characterize 
the particular model. 

For the arbitrary variables ji,d € C the scalar product of the state vectors of the model is 
given by the following expression 

{^(v±,v 2 , ■ ■ . ,v N )\^(ui,u 2 , . . -,u N )) 

N 1 1 

n^i^r n n -p±,\**q. m 

- L - L (v-jUj) M 1 - LJ - vi — vi 11 uf — ui 

j=l v J 1> N>j>k>l k J N>l>n>\ 1 n I 
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The entries of the N x N matrix Q are 

Ml -IsM U k \ M( _ X . M I U k \ I 1 



■1 • 



Using this representation and equalities (|17p that represent the steady states through the 
state vectors of the model we can calculate the projection of the state vectors (|14p and (|16p on 
the steady states. For the left steady state we have 

N 

(S N \^(u 1 ,u 2 , ■ ■ -,u N )) = lim ~\{^(vi,v 2 , . . . ,v N )\V(ui,u 2 , ■ ■ -,u N )). (19) 

{v}— >oo 

2 = 1 



This limit may be calculated with a help of the formula 
det{$(vj,u k )} _ J ^ f 1 # 



where <&(u, u) is an arbitrary differentiable for at least N times function of two variables, and 
A(-y) is the Vandermonde determinant. Taking the limit in (|19p we obtain [31j : 

N 1 
(S N \y( Ul ,u 2 ,...,u N )) = Hu 2 k ___det^ M ), 

k=l N>l>n>l U l Un 

where V^ M ^ is a N x iV matrix with the entries equal to 

vj k M) = Yx-ir ( M ) uT l ~ n \ i<3<n-i 



n=0 

M 

V (M) 
v Nk 



B-<) 

£ (-!)• ( M ) „f<»-"«). (20) 



The projection (^(ui, ii2, . . . , UAr)|SV) is given by the similar expression 

N 1 
(V( Ul ,u 2 ,...,u N )\S N ) = l[ul n -2— 2-detV-W, 

k=l N>n>l>l 1 n 

n=0 ^ ' 

-2(n-AT+l) 



n=AT 

4 Form-factors 

The calculation of the correlation function ([2]) we start with the calculation of the matrix ele- 
ments (Sn\si\^(ui, u 2 , . . . , un)} and {^f(vi,v 2 , . . . , vn)\si\Sn) which by means of (jTTJ) are ex- 
pressed through the form-factor of the projection operator, namely 

(S N \si\^(ui,u 2 , . . .,u N )) = lim {^/(vi,v 2 , . . . , ^at) |«i u 2 , . . -,u N )), 

{v}—>oo 

(^(vi,v 2 , . . .,v N )\si\S N ) = lim ($(vi,v 2 , ■ ■ • ,vn)\si\^(ui,U2, . . -,u N )). (21) 

{u}— >oo 
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By the definition (JSj) the monodromy matrix may be expressed as 

T(iA = ( B M -i(u) \ ( u Sl erf \ 

W V Cm-iH D M -i{u) ) V °f ul-u- l s x )■ 

From this formula it follows that the entries of the monodromy matrix may be expressed as 

A(u) = uA M -x{u)ax + Bm^iM^, 

B(u) = A M -x(u)<Tx + uB M _x(u) ~ u~ 1 B M ~i(u)s 1 , 

C(u) = uC M -i(u)sx + D M -i(u)crf, 

D(u) = C M -i{u)a± + uD M -x{u) - u~ 1 D M ~i(u)s 1 . 

Making use of the equalities s'j = Sj, SjO~ = vtsj = we obtain, in particular, that 

sxB(u) = (u — u~ 1 )Bm-i(u)si, C(u)s\ = us\Cm-i{u). (22) 

Substituting the explicit form of the state vectors (|14p and (|16|) into the form-factor (|21|) 
and taking into account the commutation relations (|22p we find that for the arbitrary values of 
parameters u, v the form-factor is proportional to the scalar product (1181) of the state vectors 
on a lattice with M — 1 sites: 

{$(vx,v 2 , ■ ■ . ,v N )\sx\^(ux,u 2 , ■ ■ -,u N )) 

N N N N 

= (ninc(«>ij]fl(«i)|fi> = (n| JJc , («>iII^( u *)l n ) 

j=i i=i j=i i=i 

N N N 

= n (i - u f) ("i n c'm-i(^) n ^-iwi«>- 

fe=l j=l i=l 



Taking the limit {v} — > oo we obtain 

lim 7(*(vi,u 2 ,... 

W AT 

I(i-nr; 2 ) lim (fiincM-i(^-)n^-i 

fc=l i=l f=l 



N 

(S L \sx\^(ui,u 2 , . . -,u N )) = lim ~\(^{vx,v 2 , . . . ,v N )\sx\V(ux,u 2 , . . -,u N )) 

N N N 

( Ui )\n) 



n radet^- ), (23) 



fc=l v>z>n>i u * Un 



where the entries of N x N matrix V^ M ^ are (120p with M replaced on M — 1. Respectively 
we have 



v 



(V(ux,u 2 ,...,u N )\sx\S N ) = l[ul -^—-detV^l (24) 

fc=l jV>n>«>l U l Un 

5 Bethe ansatz solution 

By following through the arguments of the algebraic Bethe ansatz method [HO [3] one proves that 
the state vectors (114p and (|16|) are the right and the left eigenstates of the transfer matrix (jlip 
with the same eigenvalues 

t(v)\^(ux,u 2 , . . .,u N )) = @ N (v, {u})\V(ux,u 2 , . . -,u N )}, 
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(*(tti,u 2 , . . . ,u n )\t(v) = {^(ui,u 2 , ■ ■ ■ ,u N )\@ N (v, {u}) (25) 
if allowed parameters u±,U2, ■ ■ ■ , un in (|2"5"j) satisfy Bethe equations [HI [TSJ HS1 El Q2J 

N 

(1 " u-T M -n 2N = ('If' 1 IK* - ("l)^ 1 ^ 2 - (26) 

i=i 

The eigenvalues &n(v,{u}) are equal to 

N 2 N 2 

®^w)=n^+(i--- 2 ) M n^- 

j=i j j=i j 



From the definition of the cyclic shift operator (|12p it follows that its eigenvalues are equal to 
N 1 

&n(1,M) = H- - 2 - (27) 

J= i 1 " u i 

From (|13p one finds the eigenenergies of the Hamiltonian ([1]) 

A' 



^ = -^e^(i, { u })^e N (v, {u})\ v=1 = -J2 -r^r 

j=i j 



The obvious solution u\ = U2 = ■ ■ ■ = un = oo provides the stationary solution with the 
eigenvalue = 0. The steady state \Sn) is the eigenvector of the transfer matrix (fTTj) t(u) 
with the eigenvalue equal to one: 

t(u)\S n ) = & N (v,{oo})\S N ) = \S N ), 

and hence it is stochastic. 

The scalar product of eigenvectors (|14[) , (|16|) is found from formula (|18p with = iij satisfying 
Bethe equations (126j) . Understanding the diagonal elements of the matrix Q in the sense of 
L'Hopitall rule one obtains the following expression for the norm of any eigenvector: 

N 2 {u!i u 2 , ...,u N ) = (y(ui,u 2 , u N )\y{ Ul ,U2, u N )) = U 2N TT det Q (28) 



with the entries of the matrix Q equal to 
N-1 + (M-N+ l)uj 2 

Qjk = : ^2 Sjk — (1 — Sjk)- 

For the special solution u\ = U2 = ■ ■ ■ = = oo the norm of eigenvectors is equal to Zjy- 

There is Zn solutions of the Bethe equations [16], and the state vectors belonging to the 
different sets of the solutions of the Bethe equations are orthogonal. The eigenvectors (|25p 
provide the resolution of the identity operator 

j _ 1^(111,112, ■ ■ . ,U N ))($(U1,U2, ■ ■ -,U N )\ 

j-f J\f 2 (ui,u 2 , ■ ■ ■ ,u N ) 

{u} 

where the summation is over all different solutions of Bethe equations (|26p . 
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6 Stationary correlation function 

For the transitionally invariant system the form-factor of the projection operator s m is expressed 
through the form- factor of the operator s\: 

F m-1„ ,1-mi 



($(Ul,U2,...,U N )\s m \S N ) = {y(ui,U 2 ,...,U N )\T m L SiT L m \S N ) 

= ll( 1_n i ) (V{ui,U2,...,u N )\si\S N ), 

3=1 



where the property ([27]) was used. 

The substitution of the resolution of the identity into (S]si\sie~^ H s m \SN) gives 

io i -\t\H ip v {S N \s 1 e~^ H \^(ui,U2,...,u N ))(^(ui,u 2 ,...,u N )\s m \S N ) 

{b N \Sie S m \£>N) = >^ TToT n 

{«} 

\t\E N (Sn\si\^(u 1 ,u 2 , . . . ,u N ))(^(u 1 ,u 2 , ■ ■ .,u N )\si\S N ) 

N / „\ m-1 



E 

W Af 2 (u 1 ,u 2 ,...,u N ) ]1 (1 -u] 



-2 



From the determinantal representations of the matrix elements (|23p , (|24p and of the norm ([28 
we finally obtain the answer for the stationary correlation function ([2]): 

^(SN\ Sl e-^s m \S N )=(^^^ 

N f„,2 i\„,2 



e -\t\E N » ( M 2 _ 1)u 2 ^ y(M-l) det y(M-l) 
^ ^ firi-.-r- 1 " detQ 



{u} j=l [I - U 

where the summation is taken over all different solutions of Bethe equations (|26|) except the 
special one. The first term on the r.h.s. is the contribution of the stationary state. 



7 Conclusive remarks 

The models connected with the 'crystal base' i?-matrix ([7]) naturally appear in the theory of the 
boxed plane partitions - three-dimensional Young diagrams placed into a box of a finite size, 
and in the theory of random walkers |32[ [5] . Some aspects of these connections were discussed 
in the papers [331 (Ml E51 EE EH [38] • 

In its turn plane partitions and random walkers are employed in analyzes of the models 
of statistical physics describing faceted crystals [39], direct percolation [3D], one-dimensional 
growth processes [HJ . It emphasizes the importance of the mentioned integrable models in the 
theory of the non-equilibrium processes. 
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